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Abstract 



We compute the partition function of the trigonometric SOS model with one 
reflecting end and domain wall type boundary conditions. We show that in 
this case, instead of a sum of determinants obtained by Rosengren for the 
SOS model on a square lattice without reflection, the partition function can be 
represented as a single Izergin determinant. This result is crucial for the study 
of the Bethe vectors of the spin chains with non-diagonal boundary terms. 
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1 Introduction 



The recent progress in the study of the open XXZ spin chains with non- 
diagonal boundary terms [H El [3l HI [5] permitted to apply the algebraic 
Bethe ansatz technique [6l [7] to this situation. It opens a possibility to 
investigate these models as well as some systems out of equilibrium such as 
ASEP. In particular there is a hope to compute the correlation functions 
using the algebraic Bethe ansatz technique [SJ El [TO] • 

In this paper we consider the trigonometric SOS model with a reflecting 
end. The reason to study such a particular case of general elliptic SOS 
model is the fact that it is related to the XXZ chain with non-diagonal 
boundary conditions by a gauge transformation [21 [3] (similar to the gauge 
transformation permitting to reduce the XYZ chain to the SOS model). Our 
main goal is to compute the partition function for this system with domain 
wall type boundary conditions. 

The domain wall boundary conditions for exactly solvable models in two- 
dimensional statistical mechanics were introduced for the first time in the 
context of the calculation of the correlation functions [TT] or, more precisely 
of computation of the scalar products and norms of the Bethe states. The 
partition function for the six-vertex model on a square N x N lattice with 
these conditions was obtained by Izergin [T2] as a determinant of an x 
matrix. This determinant representation is the corner stone for the com- 
putation of the correlation functions of quantum integrable models using 
algebraic Bethe ansatz [131 El E] . 

For the open spin chains with diagonal boundary terms the same role 
is played by the six vertex model with one reflecting end and domain wall 
type boundary conditions on the other ones. The partition functions of this 
system was studied by Tsuchiya [Mj . Once again the partition function was 
obtained as a determinant of a slightly more complicated A x A matrix. This 
representation led to a determinant representation for the scalar products 
and norms of the Bethe states [3| and finally permitted to obtain the multiple 
integral representations for the correlation functions [TO] . 

For the most general solution [TO] non-diagonal solution of the reflection 
equation [T7] for the XXZ spin chain the algebraic Bethe ansatz using the 
gauge transformation which diagonalize the boundary matrix but transforms 
the usual trigonometric R matrix (corresponding to the six-vertex model) 
into the dynamical trigonometric R matrix (solution of the dynamical Yang- 
Baxter equation [181 ITOj [20l I21j . corresponding to the trigonometric SOS 
model). 

The SOS model with domain wall boundary condition was recently con- 
sidered by different methods |221 123| . However even in the trigonometric 
case the partition function can not be expressed as a single determinant (it 
is obtained as a sum of determinants by Rosengren) and it makes the further 
steps much more complicated. 
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Here we add a reflecting end to the trigonometric SOS model. A priori 
it makes the result even more complicated. However we show in this paper 
that the corresponding partition function once again can be written as a 
single determinant. In some sense the presence of a reflecting end permits 
to avoid some difficulties of the dynamical case. We hope that starting from 
this result it will be possible to compute the scalar products and the norms 
of the Bethe vectors for the most general open spin chains. 

The paper is organized as follows. In the Section 2 we define the model 
and construct the corresponding dynamical reflection algebra. In the section 
3 following [llj we establish the properties defining the partition function 
and we show that there is a determinant solution for this functions. 



2 SOS model and dynamical reflection equation. 

The SOS model is a two dimensional statistical mechanics lattice model 
which can be defined in terms of a height function. Every square of the 
lattice is characterized by a height 9 and the its values for two adjacent 
squares differ by rj. There are 6 possible face configurations 
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and the corresponding statistical weights can be written as an R 
matrix acting in a tensor product of two two-dimensional spaces, 



R{X;t 



/RXXiX; 



V 
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i?+Z(A;0) Rt+{X;e) 

R^tlX;e) RZl{X;e) 

RZZ{X;e)/ 



If the statistical weights do not depend on the height 9 the model becomes 
equivalent to the six vertex model. 

In general the model is exactly solvable if this R matrix satisfies the 



3 



Dynamical Yang-Baxter equation (DYBE) 



RuiM -\2;9- wDRniXi - A3; ^)i?23(A2 - A3; ^ - ) 

=R23{h - As; 9)Ri3{Xi - As; - r,a|)i?i2(Ai - A2; 6) 

(2) 

The most general solution of this equation can be written in terms of 
elliptic functions but here we consider only the trigonometric solution with 
statistical weights 



( A; 0) = RZZiX-9) = sinh(A + ry) 



RXziX;e)=RzUX; 
Rt^{X;e) = R-t{Xz 



sinh A sinh(^ — rj) 

sinh^ 
sinh rj sinh(^ — A) 

sinh 9 



(3) 



We consider this model with a reflecting end, which means that each 
horizontal line makes a U-turn on the left side of the lattice. It produces 
two following configurations characterized by the weig hts }C^{X;9): 





JCX{X;9) 



/C:(A;( 



Figure 1: Boundary configuration with external height 9. 



It's important to note that such refiecting end imposes a constant ex- 
ternal height 9 for the left side of the lattice. Thus to preserve integrability 
the diagonal boundary matrix /C(A; 9) should satisfy the usual reflection 
equation 

i?i2(Ai - A2; 0)/Ci(Ai; 0)i?2i(Ai + A2; 9)1C2{X2; 9) 

=/C2(A2; ^)i?i2(Ai + A2; 9)lCi{Xi;9)R2i{Xi - A2; 9) (4) 

The diagonal solution of this equation (which corresponds to a general so- 
lution for the six- vertex ii- matrix after a gauge transformation [3j) is 

/ sinh(e+c-A) „ \ 

sinM^K+A) ^^^^ 
\ ^ sinh(C+A)/ 
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This reflecting end lead to different parametrization of the weight if they are 
in the two different half the rows. Indeed, parametrization should respect 
row and line multiplication, and also some fundamental symmetry of the R 
matrix (as the ice rule (6) , that will be presented later) . We can easily check 
that this convention lead to a well defined inhomogeneous model: 
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The domain wall boundary conditions can be easily derived from the 
corresponding conditions for the six vertex model using the one to one cor- 
respondence between the configurations of the two models in the limit of the 
height-independent weights. 
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Figure 2: Domain Wall Boundary Conditions 



One easily obtain the following boundary condition: the heights decrease 
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from left to rights on the upper boundary, the heights grow from left to 
right on the lower boundary. As left external height is fixed these two 
conditions determine completely the configuration on the right boundary 
(heights decreasing in the upward direction). 

In this paper we compute the partition function of this model. Once 
more time we would like to stress the point that while the model described 
above seems to be quite artificial it has a direct relation to the spin chain 
with non-diagonal boundary terms and this partition function is a necessary 
step to compute the correlation functions for such chains. 



2.1 The bulk weights and their symmetry 

Prior to the computation of the partition function we need to establish some 
properties of the trigonometric SOS R matrix and of the corresponding 
monodromy matrices. 

This R matrix satisfies four important properties: 



1. Ice rule 

[al + al,R,2i\;e)] = (6) 

This symmetry is responsible of the six vertex texture of the statistical 
weight: = unless a + f3 = ji + v. It is easy to see that this relation 
induce the following similar relations for the transposed R matrix 

[at-al,R'M><;e)] = (7) 

2. Unitarity 



i?i2(A; e).R2i{-X; 9) = - sinh(A - rf) sinh(A + r])Id (8) 
3. Crossing Symmetry 

The crossing relation for the dynamical i?-matrix are not as simple as 
for the vertex type R matrices, here we write it in the following compact 
form: 

- al : R%{-\ -r,-e + rjal) : <^f^^^^§^^ = ^21 (A; 6) (9) 

where we assume the following normal ordering: the in the argument of 
the R matrix (which does not commute with it) is always on the right of all 
other operators involved in the definition of R. 
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2.2 Bulk monodromy matrix and double row monodromy 
matrix 

The bulk monodromy matrix defined as 

N 

ro(A;0) = i?oi(A-ei;^-r/5^af)...i?o^(A-e^;e)= (^^l^i^J ^l^'.^Jj 

(10) 

satisfy the dynamical Yang-Baxter algebra 

N 

i?i2(Ai -\2;e-vYl erf )ri(Ai; e)T2{\2; e-iial) 



i=l 



--T2{\2; omxi-e - wDRuiXi - A2; e) (11) 



To describe a reflecting end we introduce the double row monodromy 
matrix [7] 

^(^^ «) = (c(A; S P(A;?)) = ^<'^ ""^'^^ "'^('^ 
=i?oi(A - 6; ^ - 7?^ af )...i?oiv(A - ^) 

N 

x/C(A; ^)i?7Vo(A + 0)...RioiX + 6; ^ - f^l) (12) 

As the /C matrix ([5]) solves the (ordinary) reflection equation this double-row 
monodromy matrix satisfies the following dynamical refiection equation 

N N 

Ru{\i - X2;0 -vYl <yt)Ti{\i;e)R2i{\i + As; ^ - 7? 5^ ^^f )r2(A2; 6) 

i=l i=l 

N N 

= T2{\2\e)Rx2{Xi + As; 9 - r? af )ri(Ai; ^)i?2i(Ai - As; ^ - vYl ) 

i=l i=l 

(13) 

This equation contains the commutation relations for the generators: .4(A; 6), 
B{X; 0), C(A; 6) and I'(A; 9), the only one which is important for the compu- 
tation of the partition function is the relation for the B operators 

B{Xi; e)B{X2; 9) = B{X2; 9)B{Xi;9) (14) 

It is also important to establish some symmetries of the B operators. 
Using the crossing relation (9) and the Ice Rule for the transposed R matrix 
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(7) we get: 



N 



f (A; e) =Rno{^ + ^n; e)...Rio{X + ^i; - rjj^ 



i=2 



sinh(0) 

=7(A)r-i(-A;0) (15) 



with normalization coefficients 

AT 

7(A) = (-1)^, 7(A) = (-1)"^ n + ^^-^) sinli(A + + (16) 

i=l 

It implies for the B operators: 

B{\- 9) =7(A) (tCZB{X; e)A{-\ -r];e + ri)- )C+A{X; e)B{-X - ??; - r?) 

sinh(g-7?^^^af) 

^hlh(^) ^^^^ 

And using the dynamical Yang-Baxter algebra for the bulk monodromy 
matrix, this leads to the following symmetry of the B operators: 

^ " ' ^sinh(2A)sinh(A-C + ??)sinh(A-6l-C + r?) 



(18) 



3 Partition Function 



The partition function of the SOS model introduced in the first section can 
be written in terms of the double row monodromy matrix (this construction 
is parallel to the corresponding six- vertex partition function [1^) 



i=l j=l i=l i=l j=l 

N 

={0\l{B{Xf,e)\0) (19) 



i=l 



where |0) is the state with all the spins up and |0) is the state with all 
the spins down. We will follow the standard way to compute the partition 
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function |lHll2j. first we establish a set of properties defining it in an unique 
way and then we will propose a determinant formula which satisfies all these 
conditions 

The partition function (jl9p satisfies the following properties: 



i) For each parameter Aj the normalized partition function 

^7V,27v({A}, {ei, e) = exp (^2N + 2) ^ A,^ 

X sinh(0 + C + \i) sinh(^ + A,)Zjv,2w({A}, {C}, 9) 

(20) 

is a polynomial of degree at most 2N + 2 in e^'^^ . This property follows 
immediately from the definition of the double row monodromy matrix. 

ii) For = 1 the partition function is just a sum of two terms 

sinh??sinh(6'-7?) 

A,2(A,5,6'j = —J- 

smh 

as there are only two configurations possible. 

iii) Z7v,2Af ({A}, {^}, 0) is symmetric in Aj. This property follows from the 
commutation relation (|14p for the operators B. 



iv) Zn^2n{{Mi {'^}> ^) is symmetric in ^j. This is a direct consequence of the 
Dynamical Yang-Baxter Equation. It is sufficient to insert i?j+i^i(i^j+i — 
^i\0 — 'n'^jLi+2) C21) to get the symmetry for any elementary per- 
mutation o ^j+i. 

v) Crossing symmetry. 

For any parameter Aj using the symmetry of the B operators (jlSp we 
obtain the following relation 

7 ( \ s\\sc\a\ I X ^ sinh(2(Ai ^ rj)) sinh(Ai + C) 

^iV,27V -Aj - ?7, {A}, {CI, Q) = -7 A ■ . — — ^ 

smh(2Aj) smh(Aj — C + ??) 

sinh(Ai + C + I 



. . . ^JV.2iv(A„{A},{C},{n) 

smh(Aj — d — r]) 

(22) 
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vi) Recursive relations. 

There are two points where we can easily establish recursive relations, 
fixing the configuration in the lower right or the upper right corner by 
setting Ai = ^1 or Aat = — ^i. It is easy to see that it leads to the 
following recursive relations 



Zn,2N ({A},m,W) 

N 



Ai=6 



sinhr7sinh(C — Ai 
sinh(C + Ai) 



X 



n""^(^-+^^^ sinh(g + (iV-2. + l),) 



N 



sinh(Ai - ii + ri) sinh(Ai +ii + -r]) sinh(Ai - Ci + ^) 



i=2 



X ^(Af-l),2(7V-l)({A}2...Ar, {?}2...Af, {0}) 



(23) 



Zn,2N ({A},U},W) 



sinh?7sinh(0 + C, ~ ^n) 
sinh((9 + C + Ajv) 

Ajv=-5i 



^ , sinh(g+(jV-20r/) 

n - ^^)sinh(0 + (Ar-2z + l)r7) 



X 
i=2 

X V-i),2(iV-i)({A}i...iv-i, {C}2...iv, {0}) (24) 



sinh(AAr + + f?) sinh(AAr - + ??) sinh(Ai_i + 6 + f?) 



Lemma 3.1 The set of conditions i)-vi) uniquely determine the partition 
function Zn,2n{{>^}, {C}, i^})- 

To prove this lemma it's sufficient to observe that the normalized partition 
function ()20p is a polynomial of degree at most 2N + 2 in each parameter 
e^^' defined in AN points. Indeed, due to the symmetries iii) and iv) the 
recursion relations vi) can be established for any points Aj = it^j . Due to 
the crossing symmetry v) similar recursion relations can be established in 
the points Ai = — rj. Hence we can prove by induction starting from the 
case N = 2 that the partition function is uniquely determined. □ 

It means that if we find a function satisfying the above conditions it is 
the partition function. 

Theorem 3.1 The partition function of the trigonometric SOS model with 
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reflecting end can be represented in the following form 
N 

Yl sinh(Aj + ^j) sinh(Aj - ^j) sinh(Ai + + rj) sinh(Ai - + rj) 

i,i=l 

^ 

Yl sinh(^j + ^j) sinh(^j - .^j) sinh(Aj - A^) sinh(Aj + Aj + r/) 

l<i<j<N 

(25) 

where the N x N matrix Mij can be expressed as a sum of two terms: 

_ sinh(g + C - A.) ^ sinh(C - A.) 
'""''^ - smh{e + C + A,) + sinh(C + A,) 

_ 1 1 smh((9=Fr/) 



''■^ sinh(Aj =F (,j + ^) Vs™h(Ai ± ^j) sinh ^ sinh(Aj ± + ??) 

'(26) 

To prove the theorem it's sufficient to check the properties i)-vi). The sym- 
metries follow directly from the determinant structure and the = 1 case 
is evident. The recursion relations can be easily obtained by straightfor- 
ward computation as the determinant term is reduced to a determinant of 
a (iV - 1) X (TV - 1) matrix. 

It can be convenient to express the matrix Mij in a slightly different 
form: 

^ sinh(g + C + gj) sinh(C-gj) 

sinh(0 + C + Xi)' sinh(C + AO 
^ sinh(2A0 sinh 7? 

sinh(Ai - + r/) sinh(Aj + + r/) sinh(Aj - ^j) sinh(Ai + ^j) 

In this form the result becomes very similar to the six vertex case |14] . 



Conclusion 

The main result of this paper is the fact the partition functions of the SOS 
model with domain wall boundary conditions can be in some cases (namely 
if there is a reflecting end) expressed as a single determinant. This result 
opens some new possibilities for the study of the scalar products of the 
Bethe states and finally of the correlation functions for the spin chains with 
non-diagonal boundary terms. 

There are few other interesting questions arising from our result. First 
of all it should be quite easy to generalize it to the most general elliptic 
SOS model. The result should be once again a single determinant. Also it 
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is interesting to apply it in the particular ice-like point rj = to study the 
three color model with a U-turn, following the analysis of Rosengren |22j 
and Kuperberg plj . 
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